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^BSTEACT 


Rolling pullout is one of the important, taaiical 
maneaYnns which can be used in offensive as well as 
defensive roles. The problem of designing an optimal 
ocntroller for performing this maneuver with a modern 
fighter aircraft is formulated* The effect of inertial 
cross coupling terms and the need for their inclusion in 
a mathematical model of a modern aircraft are clearly 
brought out. A five degree of ' freedom nonlinear model 
including these coupling terms and other nonlinear aero- 
dynamic terms is proposed for the maneuver problem. The 
maneuver problem is transformed to an optimal control 
problem, with the formulation of a suitable performance 
index and state model. The Gradient me cho4 of solving 
the optimal control problem and some computational 
aspects peculiar t- this problem are discussed* The 
problem is solved f^r three different flight conditi.^ns 
for a modern swept wing 'fighter aircraft and the 
resulting op'Limal controls and optimal trajectories are 
analysed. 



CHAPTER I 


IHTRODUGTIOH 

The Control problems associated vi/ith modern airc-raft 
are Qpite diversified in nature. This depends in general on 
the type of the aircraft viz., fighter, bomber,- transport etc 
and on the speed, whether subsonic, transonic or supersonic. 

In early days, autopilots were designed mainly for 
transport aircraft with a view to maintain the attitude and 
heading of the aircraft- But 'with the advent of high 
performance jet aircraft, the tactical requirements cf high 
maneuver. ability invariably ended up in a design which had 
many control problems. This is accentuated by the fact 
that stability and maneuverability are two c-^nfli-cting 
requirements in airframe design; .-ne has to be sacrificed 
for the benefit of the other. The tactical requirements 
give preference to maneuverability over stability at the 
airframe uesign stage. But stability of an aircraft is 
equally important and this is taken care of by suitable 
automatic control systems. Thus automatic control systems 
have become an integral feature of the aircraft. 

Basically the aircraft control systems are of two 
kinds ; (1) To assist the pilot in handling the aircraft 
efficiently without strain, by stabilizing any unstable 



■fflotijns. (2) To perf':nn a specified task or maneuver. 

The conventional autopilots,- pitch dampers and yaw dampers 
belong to the f-ormer classj and Instrument landing .systems: 
and Maneuver demand systems belong to the latter. The 
maneuver demand systems are those which perform a desired 
raaneuver in a specified way. In the case of transport 
aircraft, this may be an automatic landing system. In a 
fighter, this system may be designed for performing a 
complicated maneuver which the pilot may not be able to 
perform as efficiently as a control system, which can have 
decidedly better a c c u r a cy , _ j ud g em en t and speed of response. 

Traditionally, the aircraft control systems have- 
been designed using classical transfer function methods, 
where the overall reduireraents such as rise time, bandwidth 
etc. con be satisfied in ever so many ways and it was left 
te- the discretion of the designer to choose the system he 
feels best. On the other haind, application ofM.odern 
c^:ntrol techniques has distinct advantages in the design 
of a c-:-.ntrol system for performing a tactical maneuver. 

This is due to the fact that the requirements and various 
physical constraints of the maneuver can be very precisely 
specified iind the application of optimal control theory 
results in a performance which is optimum in the sense 
that the performance cannot be improved up-on any further 
without altering the problem itself. 



1 .2 SCOPE OF THE THESIS 


The thesis problem is concerned with designing an 
optimal control system for performing a very important 
tactical maneuver, namely the rolling pullout maneuver. 

The rolling pullout is one of those maneuvers 
which can be used in both offensive and defensive roles.. 

In defensive role, it may be reOiuired for evasion and 
escape from an attacking fighter aircraft at close range. 
In offensive role, this maneuver appears to be necessary 
during attacks with high closure rates,, in order to make 
aiming corrections and to achieve satisfactory breakaway 
after release of the armament. The maneuver is usually 
initiated frc.m a straight and level flight. At the end 
of the maneuver, the aircraft has rolled through about 
80° and is subjected to high normal acceleration. This 
puts the aircraft in a tight turn which is the object of 
this maneuver. The maneuver is normally performed by the 
pilot by first rolling the aircraft Quickly to one side 
through the desired angle by a step movement Of the ailero 
and then applying a nose wp elevator 'which puts the 
aircraft in a tight turn with high normal acceleration. 
Obviously, a maneuver like this, inv.,.lving the opera tiun 
of two control surfaces, can be performed in innumerable 
ways; the optimal way being the one in which the maneuver 
is performed in the shortest time with maximum normal 



4 


Tho object of this thesis is to find an optimal 
e..ntr.,l law v/hich will rOll a given aircraft through a 
desired angle in rniniiTiUra tiwe and simultaneously achieve 
niaxiraun'i normal accelercti^-n at the end of this time. This 
maneuver problem is formulated in 'this thesis as a bicriterion 
optimal control problem subjecx to (1) censtraints on the 
aircraft parameters like angle attack, side slip angle, 
rwll rate exc., (2) c -nstraints -m availoble control power 
and (3) a, specifieu attitude for xhe aircraft at the end 
of the maneuver. 

The basic Hjathematical mouel used f-^r xhe aireruft 

V 

fo-r this problem is that due tc Ehoaclsj 1 I'with a few 
simplifications. Since the model chosen is nonlinear, the 
resulting opxinial control problem cannot be solved analy- 
tically. Hence, a numerical method of solving the pr-hblem 
is indicated. Using this Mjexhod, an ezample problem is 
s.lved f^r a swept wing fighter aircraft whose stability 
derivatives are ''■aiQ\m for three flight conditijns. 
Implementation of the resulting optimal c^^ntrol law and 
cztensi-n of this method to take care of other flight 
G-uiditiOns are discussed. Sc.pe f.^r furcher research in 
this area is also iodic a ted towards the end. 

1 .3 BRIEF OUTLINE 0? THE gHlPTmS 

^ brief description of the contents of the 
vari -US c'rrjptcrs is presented below. 



Chapter II discusses the effect of inertial cress 
c-upling in nodern aircraft and emphasizes the need for 
including the terms due tu inertial coupling, in the 
eduaciwns of rco-tion of the aircraft. The effect Of 
neglecting the c jupling terms on the aircraft response 
is brought ^ut by plotting the aircralt response during 
a roll maneuver with and withc-ut these terms. 


The physical roQuirements of the rolling pullout 
maneuver and the associated control problem are stated in 
Chapter III. The chapter presents a workable model for 
the aircralt to suit the maneuver. The various assump- 
tions ■i.made ;nd the justifications in arriving at such a 
vvtrkable raOt-el are discussed. 

In Chapter IV,*^suitable performance criterion 
is formulated and the aircraft maneuver problem is 
transformed into an opximal control pr-.blem with the hehp 
•of r proper state model. The gradient method of solving 
the optimal Ci.ntrol problem is indicated. 

S -me c.,)iDputati-nal aspects peculiar to this 
problem are brougk.t out in Chapter V, The computer 
flow chart and functions 'of various subroutines are 
presented. The maneuver problem is solved for a modern 
fighter eircre.ft for three different flight conditions 
and the results are plotted. A few Conclusions based on 



these results are given. The implementstion aspects 
of the optimal control end extension of the method to 
Other flight conditions are also eiscussed. 

The thesis concludes w/ith an indication of the 
topics for further research. 



7 


CHAPTER II 

INERT IAL GOUPLIITG AMI) ASSOCIATED CORTROI, PROBLM 
IN A MODERN RICHTER AIRGRAPT 


2.1 INTROT uOT I ON 

With the increase in speed of the aircraft, especially 
after World War II, the structural and c-ui-crol pr-^blems associated 
with the aircraft have c..,.rrespundingly became mure complex. One 
such pr'.,bltrn is that due t-.^ the inertial coupling between lateral 
and longitudinal modes in a high speed fighter aircraft. Due to 
this phenomenon, applicatio-n of rolling iQ 0 ;ment to a modern 
aircraft results in some yawing and pitching moments as 'well 
which, in many cases may cause divergence in pitch or yaw 
depending on the magnitude of rolling moment. 

In this chapter, a brief review of the research done 
:n the dynamics of an incrtially coupled aircraft is given. 

The effect of neglecting the coupling terms in a maneuver 
involving rolling, are clearly br^'Ught out by comparing the 
plots of response of a modern fighter aircraft, with and without 
the coupling teriTiS taken into account. 

2.2 INERTIAL CROSS COUPLING 

'^fith increase in speed, the aircraft configuration 
changed considerably from its prewar shape, to one with short 
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slender vvings with high sweepback and with most of its mass 
c viiccntrated in the fuseloge. The shift of weight from 
wings to the fuselage resaltcd in a decrease in moment of 
inertia of the aircraft about the longitudinal axis or the 
so called .X axis, and the inertia about the lateral and 
directional axes increased. This phenomenon increased the 
coupling between lateral and long itudinal equations of 
m.. ti' n ,and con be seen by examining the basic moment 
eQuaticns of an aircraft which are, 

L = pl^ - - ly) - P'l^xz 

M = qiy - P^(Iz " ^xz 

F = fl^ - pl^2 + pq(ly - I^) + qrl^^ (2.2.i> 

vi/here L, M, F are rolling, pitching and yawing moments 
respectively; p, q, r cond I„ are angular velocities 

jL ^ Zi 

and m-'me-nts of inertia respecGively ab^ut OX, OY and OZ axes, 
is product of inertia tern. 

■ As 1-v becomes much smaller than I and I , the 
■X. j z 

m...ment of inertia difference terms (l„ - I^) and - T ) 

z X y X 

become large. Application of rolling moment to the aircraft 
results in Sume yawing mornont and th'e term pr(l -I ) may 
become large enough t.. cause an uncontrollable pitching 
raciTient. Hence this phen-'-'inenon is also known as roll 
Coupling. Moreover, 'when an aircraft rolls about an axis 
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ether than velocity vector, for example about the longitudinal 
axis with s..me angle of attack, the angle of attack after 90° 
roll is a side slip angle; after 180° of r^.ll it is a negative 
englv of attack and aftor 270° it is again a side slip angle. 
This interchange of single of attack and side slip angles causes 
the atrodynainic moments to alternate between pitching and 
yawing moments which aggravates the inertial cross coupling. 


The subject has received considerable attention in 


literature] 2-11 


Philips'; 2 Ip-nsiders a two degree of freed 


J 


cm 


aircraft rm^del including inertial coupling terms. His work 
is rasiinly concerned vvith the stability of a rolling aircraft. 


He derives expressions f^r freduencies of oscills;ti.--ns of 
rolling aircraft and vv^) in terms of the steady r^ll 

frecuency and the natural frequencies of oscillations in 
pitch and yaw of non-r-’lling aircraft and • The 

following are the results obtained by him for a rolling aircraft 
O'ith l..,wv I^. 


(1) When b:. th and greater than P^, the 

DO. ti~.ri is stable. The freciuencies of ..scillations of rolling 
aircraft under this condition are, 


w 


1 0 "nvA 


« w 

n6’ n^)- 

t J 


. r ) if ^ne ^ ^nv 


V 


(2) When either or w.^ equals and the other 

^ n0 o 

greeter than Pq, the rolling aircraft is neutr.ally stable in 
one mode . 
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(3) either of w„^ or w^,. . is less than ?. and 
the .-.-ther is greater than , the r-lling aircraft tiecoraes 
statically uns table in one mode and performs a straight 
divergence. 

(4) When both and are greater than zero but 
less than 'P^, the rolling aircraft is again stable. 

(5) When both and are much less than P^, 

V'hich means staric lo.ngitudinal and directional stabilities 
appro£:ch zer .., the two fr t.-<4,aencies of rolling aircraft 

and vV 2 b-th approo.ch . physically this means that the 
rolling aircraft can have its axis tilted frvra flight path 
and due t^ its lack of soatic stability, will continue to 
roll ab'^ut this tilted axis. 

Pr .m the above results, the following ci^nclusions 
can be dravm. 

(1) Instability may be caused by very high roll 
rates in modern aircraft with low I^. This instability 
If sts --nly as long as the aircraft rolls and would not 
therefore cause difficulty in normal flying. 

(2) The roll rate for a given aileron deflection 

romains approx im ately constant for a given True sir speed. 

But and decrease as the Indicated air soeed is 

XIB 

reduced with increase in altitude. Hence rolling frequency 

y,;-. 

may exceed one jr “other oscillation frequencies at high 


altitudes . 
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( 3 ) The ebove mentioned instability would not be 
present if pitching and yawing frequencios were equal* But 
this is an unrealiaable solution in any aircraft since the 
longitudinal stability inevitably varies with change in 

G .G .posi ti>-n. 

(4) The rate of divergence does not cause serious 
attitude change uptj 90>o angle of bank, but may cause 
serious changes in a complete 360° roll. 


Rhoads and Charaj^l ^ave carried out exhaustive- 
research on the effects of inertial cross coupling on tail 
loads during a roll pullout maneuver. The model derived by 
them is a six degree of freedom m-del which can be applied 
in general to any large disturbance maneuver. They analyse 
the r-..ll pullout maneuver tli-oroughly for the RSO aircraft. 
Also in their paper on Airplane dynamics large disturbance 
maneuvers, Rhoads and Schuler have presented a simplified 


model with five degrees of freedom to the dynamics of large 
disturbance maneuvtrs with particular emphasis on the non- 
linearities associated with high roll rates.f” 1ol 

L J 


The equations as given by them are as follows; 


* RSO A is a conventional subsonic fighter aircraft with a 
wing loading of 40 Ib/sq.ft., aspect ratio 6.4 and a 
leading edge sweep back of 9°. 
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= Z^=< - + q 

3 = p=< + YgB - r 

p = Lg(-X)e + LpP + qr (ly-I^)/!^ + (2.2.2) 

q = + M^q + pr Cl 2 “^x^^^y 

i- = Np3 + Nj.r + Np(^)p + ^ 

where the inertial cross coupling teriiis and all aajor non- 
linear terms are included. The assumptions made in arriving 
this model and the justifications in using this model for 
the rolling pullout maneuver are discussed in detail in 
Chapter III. 


Blakelock, while considering the 
cr -.ss coupling .-mits the roll transients 
rolling moment e'fuaxion and arrives at a 


effects of inertial 

r n He drops the 
I— —4 

highly simplified 


linearized four degree .jf freedom equations for an inertially 


cr^ss coupled aircraft perf 'jriaing a steady r-11. The 
equations used by him are; 


a: = q _ p p 

■J 

8 = r + ? . °< 

vJ 

( 2 - 2 . 3 ) 


where P. is the steady roll rate 


and others, 


-t 


8 Jinvestigates the 


He, like Phillips, Kclk 
stability u f the aircraft 



by applying Rough's test to the 4th degree characteristic 
e«|uation. If an aircraft is unstable during a roll, the 
stability can be iEpr^ved by increasing the "damping in yaw" 
and "damping in pitch" . This of course can be achieved by 
increasing the areas of vortical and horizontal sxabilizers 
respectively; but this calls for extensive modification to 
the aircraft and v.ill be more expensive than incorporating 
an automatic contrcl system. Mo.rewver, the increased 
stabilizer area results in increased weight (the increase 
in weight being more than that cf a control system) which 
affects the :.vcrr;ll performance of the aircraft. Blalcelock 
indicates hov/ this increase can be effected by a "^iontrol 
system inc^.rporating yaw and pitch rate feedbacks only. But 
the tactical use of this c ntrol system is 'very limited, 
as during a combat, the fighter pilot is not interested in 
rolling Continuously for more than 90° at the maximum. 

Of course this control system is useful, when ine wants to 
perform more than tvoo or three rolls, as during display 
flying . 


Thus one can see that the stability aspect of cross 
c-upling has received consiaeroble attention in literature. 
But quite apart fr^,m instability, the problem of, response 
for a stable motion may be most important for many aircraft 
during practical rolling maneuvers. In fact some aircraft 
may not be capable cf rolling fast enough (relative tw 
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their nun-rolling natural freq.uencies) to get into the 
crixical divergent region. Still they may experience 
severe and e^me times even hazardous response characteristics 
at r^ll rates bolevi? the critical value, as shovvji by 'leloh 


and 


Wils.-n 



They discuss the effect of various terms 


in the five degree freedom -oluatiuns. They also develops 


an autuiQaitiG c ntrol system for an ideal roll maneuver 


initiated by an' abrupt aileron deflection. The ideal 


maneuver described by then is one in which side slip, 
side slip rate, change in angle of attack and .angle of 
attack rate are all zero. In other words, the aircraft 
is desired to roll abjut the velocity vector. Expressions 
have been derivea for ideal yaw and pitch controls. Problems 
in fflcchanizatian ,.f this c..nxrol law have also been discussed 
Although this paper gives an elegant qualitaxive approach to 
the pr-bl-em of autumatic centr.,! during abrupt rolling 
mo-ineuvtrs, unfortunately, the control law is neither optimal 
nor has it any tactical use. 


2 .3 EPFECT OP WSG-LEGTING INERTIAL QOUPIIHG TERMS 

With a view to arrive at a simplified model for 
the thesis problem, the effect of neglecting the inertial 
c...upling terms ^n the performance of the aircraft was 
analysed for tvvo types of control inputs both involving 
roll, (i) a pure aileron roll and (ii) a roll pullout 
maneuver with saturati..n type elevator control and a step 



ailer:ri input, per both these inputs, the five degree freedom 
model (2.2.2) was analysed f.jr a swept wing fighter aircraft 
flying at 0.7 Fach at sea level. The vari:-us coefficients 
of the equations are given in Table 2.1 for the example 
aircraft for three flight conditi-ns. The following three 
cases were analysed . 

(1) The inertial c:jupling term in the pitching moment 
equation (I^. - I^)/ly was neglected. 

(2) The inertial terms in b ^th pix-ching and yawing 
moment equations viere neglected. 

(3) ill inertial terms were included* 

In all the above cases, the resulting nonlinear system 
equations were solved on a digital computer and the plots of 
=<,?-, r .11 rate and normal acceleration were obtained 
(o'igs. 2.1 t 2.4). 

Prom the plots it can be seen that 

(1) for a pure aileron roll, 

(i) Inclusion of both inertial terms gives a well 

beh.aved response exhibiting interchange of =< and p. The 
maximum roll rate is 241 deg/sec and the time taken to attain 
this roll rate is about 0.6 sec. is 1.6° and is 3° 

with a maximum normal acceleration of 3G-®. 

(ii) Inclusion Of the coupling term in yawing 

moment equation (I alone results in a divergent 

X J di 
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perf crmance . The angle of attack and si<ieslip increase 
ra nctonically in the same direction, that after 4 seconds? 
the cideslip is 23.46° and angle of attack is 19.4° with 
8 c jrresp nding n -rmal acceleration of HG®. The r ->11 rate 
olsv keeps building up to a high value of about 391 . 7 °/sec 
in 3.8 seconds. 

(iii) Omissi-.n of both the coupling terns and 

c ...xsidering a linear iriudel results in a response, v/hich 

th-n'gh bounded, is far fr-m the actual response with the 

inertial terms included. The maximuia values are = 2.3°, 

rnsjc 

8 = 3.3°, maximum r^ll rate = 260°/sec and normal 

liicl JC 

eccelerati -.n 2G®. The maximum steady roll rate is obtained 
after 1.72 sec.nds. 

A coinpariSvn of the plots leads to the conclusions. 

(1) In a rolling maneuver, neglecting the effect of 
inertial c.-upling terms gives rise tc a system v.hose both 
transient and steady state response are different from the 
actual system. 

(2) The r olling time c-)nstant is larger when the 
c.upling terms are omitted. It is nearly 3 times larger 
fur the example aircraft. 

(3) The actual side slip (with coupling terms) 
builds up first in the positive direction but then settles 
d ..wn to a negative steady state value. Neglecting the 
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inertial terms^ c coipletely distorts the response. The side- 
slip in this case never becomes negative; it settles down 
at a positive steady state value. 

(4) The actual r..ll maneuver remains a positive G 
maneuver throughout , Though there are Oscillations in G, 

it never becomes negative; but neglecting the coupling terms 
increases the amplitude of G oscillations and makes it 
negative during part of the cycle. 

(5) The interchange of angles between sideslip and 
angle of attack is not maintained when the coupling terms 
ere neglected , 

(6) The steady state value of the roll rate is an 
optimistic figure when the coupling terms are neglected. 

Thus a linear model cannot satisfactorily represent 
an aircraft performing a r>-.ll maneuver. The term 
pq in the yav/ing moment equation is seen to be 

resp-nsible for the divergent nature of the solution when 
-nly this term is included. ’ The cu-'U-pling term in the 
pitching moment equati-n seems to stabilize the divergence 
tw a new steady state. 

A comparison 'jf the plots obtained for the roll 
pullout maneuver als-o reveals that there are marked 
differences in both transient and steady state responses 



between the n'.nlinesr and linear models; and many of the 
above conclusions arrived for an aileron roll are also 
applicable for this maneuver. 

2 . 4 C QNCIUSIOWS 

Inertial cress coupling plays a very important r^le 
in a m.jdern lighter aircraft performing any maneuver 
involving rolling. Analysis of the s^^stem without inclui» 
ding the inertial coupling terms is incomplete and 
prc.'vides errone'.us response. In general, omission of 
these terms results in an underestimate of sideslip, 
angle of attack and normal acceleration and an overestimate 
-.f maximum/steady roll rate. 



TABLE 2.1 : STABILITY DERIVATIVES FOR EXAMPLE AIRGRAL^ 
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CHiiPIER III 

PhYSIOiAL DESCRIPTION Ai^ D M MHEM ATI Ci4L MODEL 
FOR THE ROLLING PULLOUT M.INEUVER 

3.1 IHIROHjcriOH 

In a matheiTiatical stud'y of any physical problein, 
it is necessary to arrive at a suitable mathematical model 
which incorporates as many physical constraints of the 
problem as is feasible, without introducing undue complexity. 
In this chapter, the physical requirements of the rolling 
pullout maneuver problem are clearly outlined. A suitable 
mathematical model is presented which will represent the 
physical system under all operating conditions encountered 
in this maneuver. Ihe various assumptions and simpllfica-' 
tions made in arriving at such a model are spelt out. 

3.2 Pi-ffol0.iL RE^U I IS 0? IHE MMEUVER PBQBLM 

Rolling pullout is one of those maneuvers which can 
be used in both offensive and defensive roles. The maneuver 
is usually performed from straight and level flight. At the 
end of this maneuver, the aircraft has, rolled through about 
80*^ and is subjected to high normal acceleration. This puts 
the aircraft in a tight turn, which is the object of the 
maneuver. The maneuver requires the operation of two 
control surfaces - the elevator and aileron, which can be 
performed in innumerable ways. 'The optimal way is defined 
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by the following two conditions ; 

( 1) The time taken to roll through a desired roll 
angle- must be a cninimum. , 

(2) At the end of this time, the aircraft must be 
subjected to maximum normal acceleration. 

These two conditions have been arrived at from tactical 
requirements of gaining combat superiority ove-r the enemy. 
The thesis problem is to determine the optimal manner in 
which the two control surfaces are to be operated to 
perform this maneuver. 

3 . 3 AGCUhATS MQDIL 

This maneuver can broadly be classified under a 
large disturbance maneuver; hence the equations of motion 
based on small disturbance theory are no longer valid for 
this maneuver. Apart from the disturbance being large, 
there are other factors like inertial cross coupling etc. 
in a modern aircraft whose effects are significant in this 
particular maneuver and which call for the inclusion of a 
few terms which are normally neglected in the equations 
of motion. 

Bhoads and 0'hara|~i conducted an extensiva 

study of this maneuver in connection with estimation of 
tail loads encountered during this maneuver. They have 
proposed a six degree of freedom model for the aircraft 
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and have derived a set of dynamic equations to represent 
the. roll pullout maneuver as completely as possible. The 
model includes many physical realities peculiar to this 
maneuver and which are normally omitted in the small 
disturbai::ice linearized theory. They are : 

(1) Changes in Mach No. are included through the 
speed equation, 

(2) The products of the disturbances are considered 
significant; thus terms like wp, vr, qr etc. are included 
in the equations. 

(3) The exact relations are used for obtaining the 
angular displacement about a set of fixed axes in terms of 
the angular velocities about the stability axes. This 
me^ans, for instance, the pitch angle ^ is not equal to 

J q dt as it is in small disturbance theory. 

(4) Variations in stability derivatives due to 
angle of attack and Mach number changes are included. 

(These terms may vary from aircraft to aircraft) Also the 
aerodynamic and inertial cross coupling terms have been 
included. 

However in their derivation the following few 
simplifying assumptions have been made to reduce the 
complexity of the equations J 

(l) The effects of non- stationary flow are omitted. 
Previous work has shown that this will not cause si.gnificant 
errorF 1 2 F 
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(2) Additional degrees of freedom due to aero- 
elastic effects are neglected. This is a valid assumption 
for any fighter aircraft which is structurally quite stiff. 

(3) Engine gyroscopic effects are not included^ as 
the moments introduced by them are normally of negligible 
magnitude. 

The increase in accuracy obtained by inclusion of these 
terms is negligible as compared to the additional efforts 
required to solve a set of more complex equations. 


3.4 A PRACTIGi^ «RKAnLE MODEL 


Though the model derived by Rhoads and O^hara 
represents the aircraft accurately during the maneuver 
considered, the solution of these equations involve 
considerable time and effort even if the available machine 
techniques are made use of. As pointed out by Rhoads 
himself pio .5 these equations can be further simplified 
without highly jeopardising the accuracy of the results 
obtained. This is required so that the equations can be 
easily simulated and solved. Here, with the help of a few 
simplifying assumptions, a five degree of freedom model is 
arrived at. 


In general, it is more convenient to express 
the equations of motion in terms of the disturbances from 
the initial conditions such that the effects of a particular 
caitrol input may be readily apparent. The maneuver 
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investigated in this thesis is assumed to he- initiated 
from a straight .and level unaccelerated flight with 
rudder locked, a rudder locked case is considered since 
the m,aneuver is normally performed with aileron and 
elevator only. Under these conditions} many of the 
initial conditions reduce to zero. Extension of the method 
to take care of other initial conditions is indicated 
towards the end of the thesis. The following simplifying 
assumptions were made in arriving at the workable model : 

(1) As the forward speed change is not appreciable 
during the maneuver, especially when the maneuver completes 
within one second, the speed equation was eliminated. This 
makes a five degree of freedom model for the aircraft; and 
eliminates the terms involving Mach number variations and 
hence greatly simplifies the equations of motion. 

(2) Sin 9 = Tan 9 = 0 and Cos 0 = 1. Even with a 
maximum value of 10 degrees for ft the maximum error would 
be less than +1 percent with this assumption. However, it 
must be emphasized that this would not necessarily be true 
if the maneuver were started from diving or climbing 
conditions. Also in this thesis, wherever possible, this 
assumption is not made use of. 

(3) 8 = v/u and < = w/u. This follows from the 
small angle assuenption. Sven otherwise^ these assumptions 
can be v«ilidated in an optimal control problem by suitable 
constraints on 8 and «: as we shall see in the next chapter. 
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(4) A few aerodynamic cross coupling terms which 
do not contribute significantly to" the p erf ormancoj have 
been omitted. 

(5) Gravitational forces are neglected as they 
are normally insignificant for the speed range considered. 

are adjusted for the usual decrease 
in aileron effectiveness at large deflections. 

(7) The equations are written with respect to 
principal axes and hence the product of inertia terms lite 
1^2 drop out automatically. - 

(8) The air speed is high enough so that the wing 
angle of attack in level flight is small. 

With the above assumptions, equations (2,2.2) were 
arrived at which are reproduced here for convenience.. 

= Z^=< - pS + q 

3 = p°< + Yp,9 - r 

p = LpC‘=<)3 + Lpp + qr (1^-1^)/!^ + C3'*4vl) 

q = + Mqq + pr ^5e 

r = + N^r + Kp(°()p + N5g_('^)6a + pq (I2.~ly)/l2. 

where the nonlinear derivatives as functions of are 
given below. These functions will vary from aircraft to 
aircraft and can be developed with the help of the 
performance curves and wind tunnel test results,. The 
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functions given below- are fee a particular swept wing 
fighter. 


L^(=<) = Lb + + 1,3^2 hl=< 

Np(=<) = Np + Np^°C + 

^6a ^ ^ ^6a2o( 

^6a-^ ''6a2 ^ a^C 


\6a.\ ) 
l5al) 


C3.4.2) 


The following equations express the orientation at any 
titn.-: of the stability axes with respect to the space axes, 
in terais of the angular velocities p, q, r:.- about the stability 
axes. These angles are marked in Fig.A-l' ih Appendix A. 

d0/dt = p + q Sin 0 Tan;© + r Gos 0 Tan 0 

d0/dt = q (ios 0 - r Sin 0 (3.4.^.) 

djf'/dt = (q Sin 0 + r Cos 0)/ Cos 9 


3.5 oOi'j 01 0 oT0i'< S 

The above equations (3.4.1) to (3. 4. 3) ) represent 
a five degree freedom model for the maneuver and ■ is a 
reasonably accurate one. The equations compare with those 
of :/lelch and Vilson also &]■ However, it must' be pointed 
out that the model is not uniformly accurate for- all flight 
conditions. For exanple, maneuvers performed at', slightly 
greater than Mach 1, will result in deceleration through 
transonic region. The aerodynamic nonlinearities so intro- 
duced may have very large effects in rolling maneuvers and 
for such cases the original six degrees of ' freedom equations 


must be used. 
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CHAPTER IV 

ftPTIMiiL GQKTROL PH)3LM 

4. 1 Hi IRQ DUCTIQH 

Although the aircraft control problems have so far 
been customarily tackled by the traditional transfer 
function methods, application of the Modern Control tech- 
niques results in several distinct advantages. Pbr a 
maneuver problem of the kind described in this thesis, it 
becomes feasible- to precisely specify the required perfor- 
mance ,and take into account the varioias physical constraints 
that are present. The resulting performance with the use of 
the control law so determined^ is optimum in the sense that 
it cannot be improved upon any further without altering the 
problem itself. Hence the maneuver problem is first trans- 
formed into an optimal control problem with the definition 
of a state model and the formulation of suitable performance 
index. The physical significance of the constraints on the 
state and the control variables is brought out clearly. The 
gradient method of solving the optimal control problem is 
shown towards the- end. 

^.2 STATE MODEL 

The system dynamic equations given in (3.4.1) can 
be transformed to the state equations by defining the 
following state and control variables. 
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(l) State Variables 



.= of attack 

=c 

in 

radians 

Xg = Side 

slip angle 

P 

in 

radians 

Xg = Bate 

of roll 

P 

in 

radians/sec. 

x^ = Bate 

of pitch 

q 

in 

radians/sec . 

Xg = Bate 

of yaw 

r 

in 

radians/ sec. 

Xg = Boll 

angle 

0 

in 

radians 

Xr, = Pitch angle 

e 

in 

radians 


( 2:|;i. Control Variables 

Since a rudder locked case is considered, the only 
two control inputs are those due to elevator and aileron. 

6e(t) = elevator deflection in radians at time t. 
6a(t) = aileron deflection in radians at time t. 
Control variables u^(t) and U2 (t) are defined such that 

^2^^^ ~ daCt)/ 

where E^a_x and are maximum possible elevator and 

aileron deflections respectively. 

With the above definitions of the state and the control 
variables the state equations can be written as follows : 
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r 


— 


{. 

^'1 


2«^1 " ^2^3 + ^4 


0 0 

^^2 


^1^3 " ^6 


0 0 

CO 


Le(Xi)x 2 + LpXg + X 4 Xgay-I^)/I^ 


^ '-^ax 

X 

4 

=u 

+ XgXgd^-I^j/Iy 

+ 

/e 0 



'1 


(x-j ) 

^5 


N^Xp + h^Xg + Np(x^)x 3 + X 3 X 4 (I^-Iy)A 2 j 


" ^ ■■ 

*6 

1 

X 3 + x^XpSin Xg + XgXyCos Xg ; 


0 0 


j 

i 

x^Gos Xg - Xg Sin Xg 


0 0 


[ 

a-„,-r- - 




(4.2.1) 

where Lp(x^) = Lp 4- U^jx^ 

= *p + V'^l ■" V® 

+ Lg^3 lUjl +(I.g^^ + Lg^s^lUjI )Xj^ 

“6a(^l)= "sa * “eaS * ^“ea=< * “sa2«'’^l 

Iht SG are the nonlinear coefficients dependent on the aii'gle 
of attack x^ and the aileron deflection u^. It may be observed 
that equation (4.2.1) is of the form 

X = fiXjU) + B(l,u) u (4.2.2) 

T 

where X = [x^ Xp ... xpj 
4 . 3 PERK)RMMCE Il^QEX 

The; problem now is to find an optimal control which 


will perform the given maneuver in the shortest possible time 
while simultaneously achieving maximum normal acceleration at 
the end of the maneuver.. Hence the problem calls for 
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optimization of two criteria simultaneously. However, this 
bi- cr itcri on problem can be reduced to an equivalent scalar 
criterion problem, by associating a suitable weighting 
factor with one of the two criteria. Thus the scalar 
performance index may be defined as 

t^ 


j" = w p-A {x(tjp + ; 


dt 


C4.3.1) 


where A-|x(t^^ = Normal acceleration at terminal time. 


The exoression for the normal acceleration at time t is 


a (t) = V(‘^ “ q + PP)/g ~ Cos 0 Cos in G units (4.3,2) 

The negative si.gn is introduced before the fmction- A since 
it is a quantity to be maximized, w is the weightage 
factor (on normal acceleration at terminal time) which 
converts the bi- criterion problem to a scalar criterion 
problem. This has been introduced, since the scalar index 
consists of two dissimilar terras (one being time and the 
other, the normal acceleration). The weightage factor 
must be chosen to give a desired relative importance 
between time' and normal acceleration at terminal time. A 
larg-'r value for W results in relatively less importance 
being assi:inGd to normal acceleratiai over the transition 
time. This is due to the fact that the normal acceleration 
is a quantity to be maximized, whereas the transition time 
is a quantity to be minimized. The effect of W on the 
optimal values of the transition time and the terminal 
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nortnal acceleration are shown for the example aircrajft in 
Chapter ¥ for a particular fli.^ht coniition. 

4.4 CON STRAI N IS 

The performance- index (4.3.1) is to be minimized 
subject to the undermentioned constraint equations and 
ineaualitic s. 

(1) D ynamic Constraints 

The system state equations (4.2.1) form a set of 
dynaaic differential constraint equations which must be 
satisfied alonq the entire optimal trajectory. These 
equations arc- to be satisfied with appropriate initial 
conditions, to be discussed later. 

(2) State Variable Constraints 

Since the problem pertains to the motion of an 
aircraft in air, the optimal trajectory is obviously 
Gontined to a particular reqioa of the state space. This 
re ;lon is defined by the constraints on the various states 
as diven below : 

(i) x^ - Anjle of Attack 

Tht^ normal acceleration of an aircraft is 
proportional to the anqlc- of attack to a first order 
approximation. Hence maximizing normal acceleration implies 
an increase in ansgle of attack with time. Thus, for the 
maneuver initiated from a straight and level flight with a 
positive or zero angle of attack, X3_ should not becorae 
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negative during the maneuver.- • Also from physical consider- 
ation} each aircraft has got a maximum permissible angle of 
attack ktiowi as the 'stalling angle', fhf upper bound on x^ 
is therefore the 'stall .angle' of the particular aircraft 
C°Siax)* bounds on x^ can be expressed by the inequality. 

This can be simplified and written as? 

1 “ 1 } — 1 0 (4,4. l) 

(ii) ^2 - Side Slip .Ingle 

The side slip during the maneuver should not 
exceed the m.aximum tolerable side slip angle (S ). This 

Hi 3.^ 

can be written as 

1 140 (4.4.2) 

(iii) x^ - Roll Rate 

Though the objective is to attain the desired roll 
angle in minimum time , the roll rate should not become too 
hi.'h to cause riisorientation to the pilot . A comfortably 
tolerable maximum roll rate may be fixed at 250 deg/sec 
^Pmax^* corresponding constraint cxi Xg is then 

I 1 - 1 0 (4.4.3) 

(iv) Xg - Yaw Rate 

The yaw rate must be of the same si.gn as roll rate, 
as otherwise there will be adverse yaw which is not. desired. 
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This is taken care of by the following inequality, 

- X5(t) XgCt) ^ 0 (4,4.4) 

(v) Xr^'- Pitch jingle 

The pitch angle must be positive through out, for 
the maiieuver considered; and must be constrained to within 
about 10 cIg ,rees (Q^ax^ that the small angle approxima- 
tion of G = Sin 9 = Tan 0 is valid. Hence 

0 x^Ct) ^ 0^^^^ 

or 1 -^(t)/6cn,ax ■ ll " 1 ^ 0 • (4.4.5) 

(vi) a (t) - Normal Acceleration 

2 -- - 

Being a positive G maneuver, the normal 
acceleration should not become negative during any part of 
the naneuver. The upper bound on a^, depends on one of the 
two following factors. 

( 1) The structural strength of the aircraft to 
withstand G. 

(2) The pilot's ability to withstand G without 
get tin,-, blacked out. 

The lower value of the above two decides the upper bound 

on a . Howfver, la modern fighter aircraft, which are 
z 

quite strong structurally, the second factor usually 
deciles the maximum normal acceleration; and for a short 
duration maneuver of this kind, an acceleration of 7 to 8 G 
(a ,^o^) can be laid as the maximum the pilot can withstand. 

2jhi a..A. ^ 
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Hence the folloi//ing constraint on Xq. 

O 


11 -140 


(4.4.6) 


(3 ) Control Variable Constraints 

The only constraint on the normalized control 
variables and .ie, that their magnitudes must not 
exceed unity. This is obvious frem the definitions of 
and Uo* Hence the constraints 

Cl 


I u^ I - 1 ^ 0 
i u^ 1 - 1 ^ 0 


(4.4.6) 

(4.4.7) 


4.6 liUTIi\L CQHDIII0H6 

The control problem is solved for the particular 
case of the maneuver initiated from a strai^'ht and level 
flight. Hence, 

= 0 (4.6.1) 

This forms the initial condition for solving the differential 
c(Xistraint equations (4.2.1) as also the initial manifold 
from where the optimal state traje-ctory starts. 

4.6 TSmib.AL TIME 

Since time during which the controls are applied 
is a quantity to be minimized, the terminal time is free. 
Also, at the terminal time, the performance index is minimum 
and its rate of change with respect to time is zero. Thus 
t^ can be found by equating dj”/dt|. to zero. 
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4 .7 lEHyili^AL MMIFOLD 

Tht perforraance index (4.3.1) is to be minimized 
subject to the constraints of Section 4.4 such that at the 
terminal time given by (4.6) above, the optimal trajectory 
lies on the folloiving terminal manifold. 

Xp ( t ^ ) = 0 

XgCtj) = 0 (4.7.1) 

XgCt^) = % 

where is the roll angle desired at the end of the 
maneuver. fhe other states are free at the terminal time 
subject to the state variable constraints which are to be 
satisfied at the terminal time as well. 

4.8 OPrihlL uOnfROL PROuLM 

The optimal control problem can now be stated 
as follows ; 

find an optimal control law u'*'(X,t) which minimizes 
the performance index 4.3.1 subject to the constraints of 
Section 4.4 and simultaneously transfers the system 
described by (4.2.1) from the initial manifold specified by 
(4.5.1) to the terminal manifold given by (4.7.1). 

As a first step towards solving the problem, the 
constraints on state and control variables are absorbed in 
the performance index itself by modifying the latter with 
suitable penalty functions! 1 3l Thus the constraints on 


Ui(t) and UgCt) ara taken care of by the addition of the 
followhag terras to the cost fu«ctlc„ within the integral 


°l,i<"i - 1) Eld - U®) 


(4.8.1) 


wher c G 


1,1 ^2,2 3-rbitrarily large positive numbers 

Call'd penalties^ and E. (l - u^) 

u.^; IS the Heavis3.de step 

function with (i - US') n.^ n 

Uj_; as Its argument such that 


0 


Ej^Cl - u^) = if (1 _ u^) ^ 

= 1 j If (1 - u£) < 0 
i = 1,;^ 


(4.8.3) 

Thr constraints on the state variables are taken 
care of by defining additional state variables- I 4 ~j. 

The inequality constraints on the state variables are 
transformed to differential constraint equations as showa 
belov/ --.ind added to the state equatiohs (4.2. .i), 


Xg - ^\<^2V1nax "D'^ “'ij 

Xg - -i| D -1^ 

^10^ VPmax^ ^ -iL D {(xs/p^^^)^ -1 






(4.8.3) 


h3= {^3V0max -D® -l} 'O' 


X = -[ (2a /a - 1 1 
13 V ^zmax 


i 

■-l} D |(S 
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.’he ; 


D. 


} = .i-l If^^. }>0 


= 0 


j 0 


(4.8.4) 


With the addition of these six equations to the original 
seven state equations} the resulting sot can be written 
in the forti 


= £(?5j m) + S(^) u) u 


(4.8.5) 


who re 


•X^ = 


[^1’ ^2 * " ^is] 


(4.8.6) 


The initial conditions of the new states (Xg to x^^) are 
zero; and on the optimal trajectory} these states are zero 
again at the terminal time. Thus the initial and terminal 
manifolds arc defined for all the 13 states. 



i\ny deviation of the terminal states from the 
terminal manifold is an error and this error must be least. 
This is achievfd by the addition of weighted error square 
terms to the performance index. Minimization of the 
performance index hence implies an automatic minimization 
of these errors. 

The perform^ance index after the inclusion of the 
weighted error square terms and the control variable 
constraints is given below 
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1 


) ^ f* “ 


2 1 


2 %,3^^3f ~ 


2 ^2,2^ 2T ■‘^2d 
^ ' 2 “ ~ 


13 


'6d^ 


r t r n 2 

+ X -^, 1 ^ ^ 


(4,8.7) 




where the subscript ’ f ’ for the state denotes the actual 
value of that state at terminal time and subscript ' d' for 
the desired value. S^ and G. . are arbitrarily large 
positive numbers. 


With the above simplifications, the optimal control 
problem reduces to that of finding an optimal control 

t) which minimizes the performance index J subject 
to the differential caistraints (4.8.6) and which 
simultaneously transfers the- system from'X = 0 to the 
desired terminal manifold. Thus the problem with the 
constraints can be identified as a Bolza type problem. 

4.9 SOLUTION BY GRADIhi^T TECaaqUE 

A closed form analytical solution to a Bolza type 
proDlem is obtainable only in the special situation when 
the constraints are linear. The differential coistraint 
equicicns in our problem being nonlinear, a closed form 
solution is not frasible and hence we take recourse to 
numerical methods. Of the various numerical methods that 
are available, the 'Gradient Method ' is chosen because it 
possesses some attractive f--atur€s i 
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Cl) Convergence is very fast in t he be ginning . 

This is 1110 rs suitable for the problein on handj whero very 
^4 gh degree of accuracy is not recjuiredj hence results in 
a saving in computation time. 

(2) It is very easy to generate an algorithm for 
a problem of this complexity. 

(3) As the nature of the solution is known before 
hand, it is easier to make a reasonably good guess of the 
initial control vector. 

A brief outline of the method as applied to the 
thesis problem is given below : 

(1) Applying Pontriagyn' s Maximum principle, a 
Hamiltonian H is defined as 

HCUjUj'Xjt) ='X^C (4.9.1) 

where L^l’^2 **’ ^3 costates and 

= 1 + i - 1) E^(l - Uj_^) (4.9.2) 

(2) The CO state equations are written by 
differentiating the Hamiltonian partially with respect to 
the state vector* 

% = - dWdH (4.9.3) 

The boundary conditions for the above costate equations at 
terminal time are given by 

ACt^) = d 


(4.9.4) 
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@ = 


_ 1 


* i - Xa^® 


■" 2 ®6,6<^>=6f ' ='6d’ ■" i 


(3) An arbitrarily small step size ’ h' is chosen 
for iiitt. gration. A trial value for the optimal control 
vector u^'^(t) is assura-d. (This is arrived from previous 
experience). The state equations (4.8.5) are integrated 
in st'-rpsj in forward time using one of the numerical 
integration techniquesj the initial conditions being 

= 0. At each step of integration, dJ/dt^ is evaluated 
and the integration is stopped when dj/dt^ = 0. The 
time tj. at which this condition is met, gives an approxi- 
mate estimate of the final time. The solution state vector 
is stored corresponding to the points of integration. 

(4) Since oi^Ct^) has been computed in step 3, 

'^^(t^) can be computed from (4.9.4). Now the costate 
equations (4.9.3) are integrated in reverse time with the 
boundary conditions at the terminal time being given 

by (4.9.4). 

(6) The components of the gradient vector 
and ^lVc^2 computed from 

2 

AH/iu = M.. E„ + G, ^u_B(u.'' - l) 

on/ou^ - 4‘-^5e%aX 1, 1 1 1 


(4.9,5) 
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d H/o'Ug 


+\CS6a " 8N5a2 '"s' "<®6a« 

+ Gg^o^gSCUg - 1) 


* ®^8a2« l“i?t>*i3‘W 
+ 2H5^3« I«s 1 )xi'3Va^ 

(4.9.6) 


The iradie-xit vector dH/du will be null only when the 
conditions for a minimura are satisfied. The change in the 
performance index at N^h iteration to a first order 
approximation is 

= LoH/cu J Au (4.9.7) 

where A ue is the change in the control vector at the N"*^^ 
iteration. If we wish to make the greatest or steepest 
change .4' in J , this can be effected by making A directed 
opposite to the gradient so that the dot product of 
and is maximum. 

(i-t;-) = - k (lb‘H/o’U^J 

or a/ = - K LaH/du“3^ D'H/a/3 (4.9.S) 

where. K is an arbitrary positive constant less than unity. 
Usually K is chosen to give a lO percent change in the value 
of the performance index. 

For a lO percent change in J, K = o'H/v3'^>c’H/c'U > 

where < . > denoto-s scalar product. 

(6) The new trial value of the control for the N***!^^ 
iteration, is given by 

= u'^ +Au^'^' 


(4.9.9) 
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(7) With this new control, are 

compute, d as before and the iterative process repeated until 
the change in the cost function (.^J) between two successive 
iterations is below some specified small value. 

(8) Under J minimum condition indicated by negligible 

change in J, the terminal manifold equations will also be 

satisfied to a reasonable degree of accuracy. The accuracy 

is increased by increasing the appropriate penalty S. . and 

1 , 1 

repeating steps l to 7 once again with the optimum control 
obtained in the first stage. 

m example problem is solved in Chapter V for a 
modern fighter aircraft. 

4 . 10 CO NCLUSIONS 

The maneuver problem is formulated as an optimal 
control problem incorporating all the physical constraints. 
The problem is solved for a swept wing fighter aircraft; and 
some canputational aspects of solving the problem by Gradient 
method are discussed in the next chapter. 
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GRATER V 

SOKE GOMPUlA'nONiVL ASPECIS AND RESULTS 

1 li-iniODUCXIOH 

The special featur'-s of the gradient^ method as 
applied to the optimal maneuver problem are discussed here 
in greater detail, a block diagram of the computer flow 
chart and the functions of the various subroutines are 
given in Figure 5.1. An example problem is solved for 
different flight conditions and different values of W and 
the optimal control and trajectories are plotted and 
thoroughly discussed. Extension to other initial conditions 
and implementation aspect of the optimal control are also 
p rt S'V nt td . 

5.2 GERTiaN GOMPU f AXIQN.iL -ASPECTS 

The following aspects of computation need special 
considera ti on. 

( 1) Initial choice of control 

Since this is a nonlinear problem solved by the 
gradient method, the control vector u will be obtained as 
a function of time. Later on, after the optimal controls 
have been computed for different flight conditions, one may 
attempt to solve the so called vSynthe sis problem by 
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expressing optimal control as a function of the state and 
time* Hencoj to start with^ the control vector wg^s arrived 
at for the example problem by finding the response of the 
aircT'ift to three different types of probable control inputs 
in performing this maneuver : 

(i) A negative step aileron -and a step elevator 
applied simultaneously. 

(ii) A negative step aileron and a ramp type elevator 
saturating at the end of l second. 

(iii) A negative step aileron; and a step elevator' 
applied 0.4 second after aileron, is applied. 

Froui the roll rate, roll angle, side slip and normal 
acceleration responses for the above three eases, it was 
decided to choose, a negative step aileron for 0,4 second 
followed hy a positive going ramp type aileron saturating 
at th> '■nd of 0.8 second; and a ramp elevator saturating 
in 0.1 second- This type of control iaput is expected to 
provide a response in between those of (i) and (ii) of the 
trial inputs and ivill he close to the desired response.- 

(2) Inte grating technique 

The system equations were integrated by Adames 
fourth order Predictor corrector method. The solution at 
the first 4 points was obtained by fourth order Runge Kutta 
method with Gill's variation. The Predictor corrector, 
though very accurate, was found to he highly sensitive to 
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the step of integration for this problem. A larg-er 

step size and high accuracy resulted in oscillatiois during 
integration. This was predominant when the costate equations 
were integrated with large values for initial conditions. 

Since the basic problem is to find the iiitegrated value at 
each in te -ration point and not the final value alone, 
automatic reduction of step size as the integration proceeds, 
cannot be resorted to. On the other hand reducing the step 
size to arbitrarily small value calls for increase in computer 
memory since the solution at the integration points are to be 
stored for calculating the gradient, performance index etc. 
However, this difficulty was overcome by a few trial runs 
v/ith various step sizes and ultimately a step size of 0.0025 
sec was decided upon with 1 percent accuracy for the 
exuTiple problem. 

(3) Ter minal time 

Since the terminal time is not fixed, the forward 
iiitegration of state equations was stopped when the total 
time derivative of the performance index (dJ/dt^.) passed 
through zero. This is only a necessary condition. The 
sufficient condition, namely the second derivative must be 
positive Was not tested duo to complexity of the expressions 
involved and consequent increase in computation time. Phis 
gave rise to stopping the integration at incorrect times. 

(The inte yratioi .1 st:..-ps when dJ/dt^. is zero for the first 
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time as the integration proceeds; but this need not 
xiecc-ssarily be the final time at which many other 
conditions are to be satisfied.) These incorrect points 
were to be eliminated by examining the terminal states at 
these points which will be far away from the desired 
terminal manifold; and also the controls will be operating 
outside the boundary. In the thesis, the algorithm has 
been developed in such a way that the test for dJ/dtjf to 
be zero is done only when the states are in the viscinity 
of the terminal manifold. This avoids testing for 
dj/dt|. = 0 find d J/dt^ ^ 0 at all integration points and 
hence saving considerable computer time. 

(4) otep size down the gradient 

The gradient method as applied to this problem is 
highly sexisitive to t he step size K which decides the change 
in control required for the next iteration. Too large a' 
step size may result in an increase in performance indvC-x 
iustead of a decrease. Too small a step size on the other 
hand may increase the number of iterations before satis- . " 
factory convergence is obtained. In the example problem, 
the step size was calculated for 10 percent change in the 
value of the- performance index; the control vector was 
modified accordingly and the system equations were integrated 
and the ri< w value of the p-‘rformance index was evaluated. 

If this luw value is less than the previous one, the process 
is continued, till such time the vaiuis of the index starts 
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incrGQ-sln- 5" this stage the step siz€ is reduced to 
1 / lO of its calculated value (for 10 percent change) and 
the orocess repeated. 

(5) Penalties 

The constraints on control variables were very 
rarely violated unless the region of search for dJ/dt^- = '0 
was incorrect. The constraints on state variables and 
terminal manifold wore usually violated. This violation 
was gradually removed by increasing the penalties on those 
variables in stages. One stage of iterations was carried out 
with one set of penalties till there was insignificant 
chaixge in the value of the performance index between suc- 
cessive i t'-'-rations. Then the penalties increased and the 
whole iterative cycle started once again with the controls 
obtaii'ied from the- previous stage. Thus the optimization 
has to be performed in stages; the penalties being increased 
at the beginning of each stage. Ultimately the penalties on 

Q 

a few states were increased to 10 for the example problem* 

(6) I nitial conditions 

The initial conditions for solving the state 
equations remain constant throughout the process of optimi- 
zation. But the initial conditions for the co state equations 
is proportional to the error in the violation of the state 
terminal manifold ^and are also proportional to the penalties 
on the states. Hence either large violation or large penalty 
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r 0 suits in iB-r.-^? vdlues for the initial conditicos* ThRse 
large initial conditions resulted in instability duping 
integration by Predictor corrector method. Hence it was 
absolutely necessary that as accurate a result as was 
possible be obtained in any stage with one set of penalties 
before the penalties were increased for the next stage. 

5.3 ELpy ouahj;' 

A block diagram of the computer flow chart is given 
in Figure 5.1. The various subroutines used are ; 

(1) CO.j TfiL 

The initial choice of the control vector is either 
read through or generated in this subroutine. Various 
tolerances which require frequent change are also fed 
through this subroutine. 

(3) oOLn 

This subroutine integrates either state or costate 
equ-itions (depending on the index value given by the main 
program) by calling subroutines PREDTX or PBEDTL respectively, 
and stores the integrated result at each integration point, 
l^fhile solving state equations, change of the sign of dJ/dt^ 
is searciied in the proper zone and integration is stopped 
when there is a change in sign. 

(3) PKEbfX and PRBf il 

These are the subroutines where the state and the 
costate equations are written. 
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(4) DOrp x^D 

LhjL s subrouting' finds thg dot product of two 
vfctors. This is used to find the dot product of the 
gr-adiont vector with itself, for estimating the step size 
K down the gradient. 

(5) iiAiv 

The subroutine finds the maximum absolute value of 
a given array of real nmibers. This is required to norma- 
lize. the optimal state solution at the solution points 
be fore obtaining a plot. 

(6) PLO T 

This subroutine provides a normalized plot of the 
optifnal solution. The plot is not very accurate since the 
normalized values are truncated before plotting; but gives 
a r‘ a eon ably goo (3 ide.a about the solution. 

5.4 bl'^iPLE 

Thr rolling pullout maneuver problem is solved 
for a swept wing fighter aircraft with the following leading 
particulars : 

W./S- = 65 lb/ sq.ft. 

Aspt ct ratio = 3.5 
oweep back = 4o® 

The stability derivatives for this aircraft were known for 
three different flight conditicos. These derivatives were 
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give?!! in rnblt- 2» 1, CIie optimal control nnd. th.6 optimal 
response were found for three different flight conditions 
x-/ith value of V/ = 1. 

(A) 0.9 Mach at 20000 ft. 

(B) 0.7 Mach at sea level 

(C) 0.7 Mach at 30000 ft. 

For flight condition the response and control were 
found for Vj = 1000 and W = 0.0001 ’.alsoj to study the 
effect of the weightage factor ¥ on the terminal time and 
normal acceleration at terminal time. The results are 
plo tted. ( Figure s 5 .2 5.6). 

^ • 5 DISCUSSIO N OF THB fl£:DULfS 

Though it is impossible to generalize based on the 
results for only three flight conditions, some broad 
conclusions can b' drawn and the general trend of the 
solution discussed. With these in view, a discussion of 
the optimal control and the resulting optimal trajectory 
and the conclusions drawn thereupon may be listed as 
follows : 

( l) Iran sit ion time 

The minimum time reciuired to perform this maneuver 
is given in fable 5.1 for the three flight conditions. It 
is sfcn that for the same Mach number, the maneuver performed 
at a higher attitude (hence reduced true air speed) requires 
longer time. This is suggestive heuristically as well. 
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(T) Normal acceleration at termina l time 

Normal acceleration decreases with increase 
ill altitude. This is so because of the decrease in control 
effectiveness at high altitudes. For the same Mach number, 
the pilot flying low caii maneuver his aircraft faster and 
c?an apply a higher value of G than while flying at higher 
altitude s. 

^haturc of the optimal control 

oincc the' maneuver calls for rolling fast to one 
side, a step aileron is suggested. Also since it is a 
stop tJ stop m-ineuver, (zero roll rate at the beginning 
and r-nd of the ■.umf'-uver) at some sta5e the aileron must be 
reversfd and a negative aileron applied so that the roll 
rate bteorn: s zero when the aircraft has rolled through the 
d- sired roll angle. This arguiient concurs with the type 
of optimal coutrol obtauied. As for the elevator, since 
the normal acceleration is required to be maximum only at 
the end, maximum elevator deflection must be present only 
towards the end. This suggests a saturation type elevator 
input. A study of the results shows that the nature of the 
optunal control obtained coincides very well with the 
nature of the control suggested heuristically. 

( 4 ) Co nstraint s 

It is seen that the constraint on P- is violated 
to a certain extent towards the end of the transition time. 



This violation can be reduced by a further increase in 
the penalty on that state; but since the \aolation was 
iiOt very severe, this process was not carried on any 
further. Also the constraints on most of the states were 
fixed by pilot's opinion rating. Tnis being highly sub- 
jective; violation of the constraints by a small percentage 
need not be viewed upon seriously. 

‘■aa- ximufii roll rate 

A comparison of the roll rate plots for flight 
condition B without elevator (Pigure 9.1) and the optimal 
roll rate shjws that the maximum roll rate in the optimal 
case is less thari that obtained for pure aileron roll, and 
the time taken to roll through 80° is more in the optimal 
case. This seeias to question the optimality of the results 
obtained. But it must be borne in mind that the optimal 
problem was solved for a stop to stop maneuver, which calls 
for zero roll r ite at the terminal time. In a pure aileron 
roll, the roll rate was approximately 80 percent of the 
maximum roll rate when 80° roll was executed. 

(6) E ffect of wei^htage factor W 

It was mentioned in Chapter IV that the weightag® 
factor W determines the relative importance in optimization 
between the two queurtities (l) the transition time and 
(9) the normal acceleration. Analysis of the results 
reveals that for a given flight cendition, the effect of ¥ 



57 

is not very significant. For a change in ¥ from 0.0001 to 
1000, the corresponding change in transition time is only 
.005 secaids and the chaiige in normal acceleration is 0.02G. 

JAflOrj OF THE OP JIM COiJfRQL LAW 

Fortunately the optimal control law for this 
problem is not highly complex. It can be easily imple- 
mented by suitable hardware. But before implementing, the 
problem must be solved for various combixiations of Mach 
nuiViber and altitude and also for various initial conditions, 
fhe results obtained are to ue compared and the optimal 
control must be expressed as a set of empirical relations 
in Mach number, altitude and initial conditions. These 
empirical relations c:an be transformed into a hardware 

realization using suitable logical circuits, amplifiers, 
actuators etc. 

The maneuver is initiated by the pilot by an abrupt 
movement of the stick to one side which actuates a limit 
switch. Thereupon^ depending ai the conditions of flight, 
suitable control is applied by the logical circuits* The 
pilot takes over from the control system when the aircraft 
rolls through the desired roll angle. 

6.7 a XJEi^olOrj fO OTHBR IiaflAL OOKDiriQMS 

The maneuver problem was solved starting from an 
initial condition corresponding to straight and level 
unaccelerated flight. Although, the scope of this thesis 
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lOGs not cover initial conditions other than straight and 
level fli:;ht, for other initial conditions, the original 

can be used 'ind the simplifying assumptions are to be 
made based on the characteristics of the initial conditions. 
For exa.iple if the maneuver is performed from a true banked 
turn iii a horizontal plane, the initial side slip velocity 
is 2 t.r.o and the initial bank angle is a finite value. There 
will be finite initial pitching velocity which is balanced 
by the elevator and a constant yaw rate. Also all the three 
controls arc^ deflected ixiitially in a true banked turn and 
the rudde-r is assumed to be locked at this deflected 
position. These are soire of the aspects to be borne in 
mind and the equations of motion should he modified 
accordingly. 

5.8 G Qi'j CL u ol 0 dS 

Certain computational aspects peculiar to this 
problem are discussed and the computer flow chart and 
functions of various subroutines are given. An example 
problem in roll pullout maneuver is solved for a modern 
fighter aircraft whose stability derivatives were known 
for three different flight conditions. The resulting 
optimal coiitrol and trajectories are plotted. Conclusions 
based on the results obtained and implementation aspect of 
the ccontrol law are also presented. 


six degree of freedom equations of Bhoads and 0 


hara 
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TABLE 5.1 : RESULTS OP ?''1NEUVER PROBLEM 


ELIG^'IT 

cord:' . 

MACH 

x^TUrSER 

ALTI I'JDE 
PT. 

W 

TRANSITION 

TIMEggQ ^ 

NORMAL 
ACCLN Q 

A 

0.9 

20,000 

0.0001 

0.6950 

4.9106 

A 

U.9 

20,000 

1 

0.6925 

4.9074 

A 

0.9 

20,000 

1000 

0.6900 

4.8906 

B ^ 

0.7 

SEA LEVEL 

1 

0.6800 j 

8.2592 

! G 

j 

0.7 

1 

30,000 

1 

0.9600 I 

1 ' - J 

1 .8490 


TABLE 2 EPEEOT 01 - H'- Oil OPTIMAL GO^srPOL 

-. ILIGHI GOIMDITION 'A' 
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CHAPTER VI 
CONaUSIOHS 

The problecQ of designing an optimal controller for 
performing the rolling pullout maneuver in a modern fighter 
aircraft was formulated. The effects of inertial cross 
coupling were brought out and the need for a nonlinear model 
including inertial cross coupling terms was emphasized. With 
this tflodel and a suitable performance index, the maneuver 
problem was converted to an optimal control problem. A 
numerical method of solving this problem was presented. An 
example problem was solved for a modern fighter aircraft 
for three different flight conditions. The implementation 
aspects of the optimal control law and extension to various 
initial conditions were also given. 

Scope for further research 

Due to nonavailability of sufficient data, it was 
not possible to solve the problem for more than three 
flight conditions. With the resxilts obtained for these 
three flight conditions, the nature of the control law is 
obvious; but with these results alone it is not possible to 
define the switching points of the control as a function of 
the Mach number and altitude. If sufficient data for 
evaluating the stability derivative': for all flight 
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conditions were available, the follo'^dng topics merit 
further study : 

(1) lo find optimal control for various flight 
conditions covering the entire flight regime of the 
aircraft and try to arrive at an empirical relation for 
the control in terms of air speed and altitude. 

(2) To study the effects of various initial 
caiditions on the results obtaixied in (1). 

(3) fo design suitable logical circuits to take care 
of altitude, air speed and ixaitial conditions. 



69 


APPEIIDIZ A 

aQHVMTIOUS im PSPIMTIQIS 

The sign conventions and definitions of certain 
terms used in this thesis are given below. 

NAPA STf^BILI'TY AXES : To specify the position of a moving 
aircraft, it is necessary t: adopt an axes system. The 
axes system adopted in this thesis is according to MCA. 
Convention and is kncjwn as the ’’Stability Axes System". 
These ere three mutually perpendicular axes with their 
origin coinciding with the centre of gravity of the 
aircraft. The axes system is fixed to the aircraft and 
moves with the aircraft. OX is positive forward and lies 
in the plane of symmetry of ’uhe aircraft with OZ alsa in 
the same plane with positive downward. OY is perpendicular 
to' this plane of symnie'cry and is positive through the 
right wing. The axes system is shown clearly in 
Figure A-1 * For a given initial flight condition^ OX is 
fixed parallel to the relative wind vector at the start 
of the maneuver. 

SPACE AX.ES I This is a set of mutually perpendicular axes 
(0Z^_j, OY^, OZ^) with OZ^ and OY. axes in the horizontal 
plane and the origin fixed at a particular point in the 


space . 
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DISPLACMBUT .^yrGLES : 

( 1 ) Pitch angle 6 : It is the angle made by the 
X stability axis with its projection on the horizontal 
plane, 

(2) Boll angle 0 : It is the angle through which 
the stability axes system must be rotated about the 

X stability axis to make the T stability axis parallel 
to horizontal. 

(3) Yaw angle V ; It is the angle between the 
projectijn of X stability axis on the horizontal plane 
and the X^ space axis. 

SOME SICT gOYVMTIOKfS : 

(^) A linear displacement along a positive 
reference axis is considered to be positive, 

(2) An angular displacement or movement is 
Considered positive if clockwise when viewed from the 
origin along a positive reference axis. 

(3) Both linear and angular acceleration and 
velocities are G-;nsidered positive in the same sense as 
their corresponding displacements. 

(4) Contr-^1 surface deflections are considered 


in the same sense as angular displacements. Ailerons 
move in opp;3site directions; hence the motion of the 
right ailer'On has been chosen to determine the sign of 



■the aileron angle. Ther e-lore, right aileron duwn and 
left ailerm up are c-nsidered p-sitive and cause a 
negative roll according t-o the sign conventions for 
angular displacements. 

REM II PIT MOVING AM) FIXED AXES SYSTEMS 

Since the stability axes move in space,, to 
enable proper definition of angular velocities and dis~ 
placements of the aircraft, the space axes are to be 
rotated in a- particular sequence to orient them with the 
stability axes. The se<^’uencG is not particularly 
imp>-.rtant as long as the same sequence is used through 
out. In this thesis, yav? pitch roll sequence is used for 
orientation of the space axes with the instantaneous 
position Of stability axes. The orient a tijn is d^.^ne in 
the following manner : 

(1) The "space axes" (OX,.^, OY , , OZ^) are defined 
at the uri.gin -f the "stability axes system" as shown in 
Figure A-1 . 

(2) The space axes system is rotated about the 

axis through an angle -y (yaw) to the orientati,^n 

^2'^2^2' ^ 0^0 rotated through an angle'Y*. )• 

* 

The angular velocity ab’^ut Z,^ isY* 

(3) Prom the new orientati-'.n space 

axes system is r.jtated about the ^2 axis (pitch) 
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thruugh an angle 6 to the new position X^Y^Z^. The 

♦ 

angular relocity about axis is 0. 

(4) Prom the position assumed in (3) above, 
the axes system is rotated about X^^ axis (roll) to the 
final orientation ]CYZ. The d.isplacement angle is and 
the angular velocity about the X axis is jZJ, 

0, e, ^ IF TERMS OP p, 1, r ; 

let p,l,r be the t.:tal angular velocities of the 

aircraft (thus of the stability axes system) about the 

* • 

X,y and Z stability axes respectively. Resolving jZl, © 
and al .ng the stability axes, we get 

^ = p 4 - q Sin ^ fan © + r Cos ^ Tan © (A.1 ) 

e = q Cos - r Sin ^ (A. 2) 

■^ = (q Sin jZl 4 r C ^s j2()/Cos © 

0^, 6 and 'y can be obtained by integrating the ab^ve 
equations. 
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